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Abstract. Neutrino-induced pion production on nuclear targets is the major inelastic channel in all present-day neutrino- 
oscillation experiments. It has to be understood quantitatively in order to be able to reconstruct the neutrino-energy at 
experiments such as MiniBooNE or K2K and T2K. We report here results of cross section calculations for both this channel 
and for quasielastic scattering within the semiclassical GiBUU method. This methods contains scattering, both elastic and 
inelastic, absorption and side-feeding of channels all in a unitary, common theoretical framework and code. We find that 
charged current quasielastic scattering (CCQE) and 1 Tt production are closely entangled in actual experiments, due to final 
state interactions of the scattered nucleons on one hand and of the A resonances and pions, on the other hand. We discuss 
the uncertainties in the elementary pion production cross sections from ANL and BNL. We find the surprising result that 
the recent In production cross section data from MiniBooNE are well described by calculations without any FSI. For higher 
energies we study the validity of the Bloom-Gilman quark-hadron duality for both electron- and neutrino-induced reactions. 
While this duality holds quite well for nucleon targets, for nuclear targets the average resonance contributions to the structure 
function F2 are always lower than the DIS values. This result indicates a significant impact of nuclear effects on observables, 
reducing the cross section and structure functions by at least 30-40% and changing the form of various distributions. 
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INTRODUCTION 

Neutrino oscillations experiments search for a distortion 
in the neutrino flux at the detector positioned far away 
from the source. By comparing both near and far neu- 
trino energy spectra, one gains information about the os- 
cillation probability and with that about mixing angles 
and mass squared differences. However, the neutrino en- 
ergy, that enters critically into the oscillation probabil- 
ity, is not directly measurable but has to be reconstructed 
from the reaction products. A proper understanding of 
neutrino-nucleus interactions is, therefore, essential for 
the interpretation of current neutrino oscillations exper- 
iments. Present disappearance experiments use the 
CCQE reaction both as signal event and to reconstruct 
the neutrino energy from the outgoing muon. CCQE is 
defined as £^p on a single nucleon; in the nucleus, 

CCQE is masked by final-state interactions (FSI). Thus, 
the correct identification of CCQE events is directly re- 
lated to the question of how FSI influence the event se- 
lection. The main background to CCQE is CC 1 pro- 
duction. If the pion is absorbed in the nucleus and/or not 
seen in the detector, these events can be misidentified as 
CCQE. Consequently, a proper understanding of CCQE 
and CC 1 is essential for the reconstruction of the neu- 
trino energy. 

The main task in a Vg appearance experiment like 
MiniBooNE is to detect electron neutrinos in a (almost) 
pure beam. The signal event, the Vg CCQE inter- 
action, is dominated by background. A major problem 



comes from misidentified events, mainly because of the 
fact, that the MiniBooNE detector cannot distinguish be- 
tween a photon and an electron. Thus, induced neutral 
current n'^ production, where the tt" decays into two 7s, 
is the major source of background when one of the pho- 
tons is not seen or both Cherenkov rings overlap. 

As all of the present oscillations experiments use nu- 
clear targets, it is mandatory to consider FSI, i.e., pion 
rescattering, with and without charge exchange, and ab- 
sorption in the nuclear medium. A realistic treatment of 
the final state interactions (FSI) can be achieved in the 
framework of a coupled-channel transport theory — the 
GiBUU model. 

After a brief review of our model, we first discuss 
the impact of pion production on CCQE measurements. 
We further investigate the influence of nuclear effects on 
CCl7i;+ and NCItt" cross sections, and, where possible, 
we confront our model to recent data measured at Mini- 
BooNE. In the second part we extend the description to 
higher energies using duality arguments. 



GIBUU MODEL 

We treat neutrino-nucleus scattering as a two-step pro- 
cess. In the initial-state step, the neutrinos interact with 
nucleons embedded in the nuclear medium. In the final- 
state step, the outgoing particles of the initial reaction are 
propagated through the nucleus using a hadronic trans- 
port approach. 
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FIGURE 1. Total p — > Tt'^p cross section as a function 
of the neutrino energy compared to tire pion production data of 
of ANL (Refs. |5] (•), |6] (■)) and BNL (|7] (x)). Tlie solid 
line has been obtained with a form factor fitted to the ANL data, 
the dashed one is fitted to the BNL data. 



In the energy region relevant for MiniBooNE, Sci- 
BooNE and K2K, the elementary vA^ reaction is dom- 
inated by two processes: quasielastic scattering and the 
excitation of the A resonance (P33(1232)). Additionally, 
our model includes the excitation of 12 A^* and A reso- 
nances with invariant masses less than 2 GeV and also a 
non-resonant single-pion background. Details are given 
inlEh. 

The excitation of the resonances (R) is described 
within the isobar model with the help of the nucleon, 
N — A and N — N* form factors. Vector form factors are 
derived from MAID helicity amplitudes 121, |3i |4j] ex- 
tracted from electron scattering experiments. 

Experimental information on the N-R axial form fac- 
tors is very limited. Goldberger-Treiman relations have 
been derived for the axial couplings, but there is no in- 
formation about the dependence. We apply PC AC and 
pion pole dominance to derive one of the axial couplings 
for each resonance and to relate it to the pseudoscalar 
coupling. The dependence of the A axial form fac- 
tors are fitted to either ANL or BNL bubble chamber 
neutrino-scattering da/dQ^ data for the v^p ll^K^p 
reaction. Fig. [T] shows the integrated cross section to- 
gether with the data. We note already here, that the solid 
curve fits the ANL data, while the dashed curve fits the 
BNL data. Thus, the latter would obviously lead to higher 
pion production cross section also on the nucleus. 

The single-TT non-resonant background cross section 
C7bg includes vector, axial and also interference contri- 
butions 
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The vector part is fully determined by electron scattering 
data, as described in IH. The axial and the interference 
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FIGURE 2. Total CC pion production cross sections for the 
mixed isospin channels as a function of the neutrino energy 
compared to the pion production data of of ANL (Refs. |5] 
(•) and 16] (■)) and BNL (|7] (x)). The solid lines denote the 
our full result including the non-resonant background. Further- 
more, we show the results for pion production only through the 
excitation and the subsequent decay of all resonances (dashed 
lines) or through the A alone (dash-dotted lines). Note that the 
ANL based fit for the A axial form factor has been used here. 



term collected under the label "non-V" are only present 
in neutrino scattering and are fitted to the available neu- 
trino data for both, v^n jJ.^n'^n and v^n jj.^n'^p. 
The final results are shown in Fig.|2] 

The neutrino-nucleon cross sections are modified in 
the nuclear medium. Bound nucleons are treated within 
a local Thomas-Fermi approximation. They are exposed 
to a mean-field potential which depends on density 
and momentum. We account for this by evaluating the 
above cross sections with full in-medium kinematics, i.e., 
hadronic tensor and phase-space factors are evaluated 
with in-medium four- vectors. We also take Pauli block- 
ing and collisional broadening of the outgoing hadrons 
into account. Our model for neutrino-(bound)nucleon 
scattering is described in detail in 1 1]. 

After the initial neutrino-nucleon interaction, the pro- 
duced particles propagate through and out of the nucleus. 
During propagation they undergo FSI which are sim- 
ulated with the coupled-channel semi-classical GiBUU 
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FIGURE 3. Ratio of the CCQE-like to the true CCQE cross section as a function of the proton (pion) momentum detection 
threshold for CC on ^^C at Ey = I GeV. The sohd lines are obtained using the tracking detector identification, while the dashed 
lines are for Cherenkov detectors. 



transport model (for details, see Ref. f§\ and refer- 
ences given there). It is based on the BUU equation 
which describes the space-time evolution of a many- 
particle system in a mean-field potential including a col- 
lision term. Nucleons and resonances acquire medium- 
modified spectral functions and are propagated off-shell. 
Herby we ensure, that vacuum spectral functions are re- 
covered after leaving the nucleus. The collision term ac- 
counts for changes (gain and loss) in the phase-space 
density due to elastic and inelastic collisions between the 
particles, and also to particle decays into other hadrons. 
Baryon-meson two-body interactions (e.g., 7zN 7zN) 
are described by resonance contributions and a small 
non-resonant background term; baryon-baryon cross sec- 
tions (e.g., A^A^ NN, RN NN, RN R'N, NN 
kNN) are either fitted to data or calculated, e.g., in pion 
exchange models. The three-body channels nNN NN 
and AA^A^ NNN are also included. The BUU equa- 
tions for all particle species are thus coupled through 
the collision term and also through the potentials. Such a 
coupled-channel treatment is required to account for side 
feeding into different channels. 

The treatment of pion final state interactions in 
GiBUU has been widely tested both with nA and 7A 
nX data |9, 10, 11]. The latter reaction is quite similar to 
the vA nX reaction in that the incoming particle inter- 
acts with all target nucleons and the vector couplings are 
the same. 

CCQE AND CCl;r+ ENTANGLEMENT 

One challenge is to identify true CCQE events in the de- 
tector, i.e., muons originating from an initial QE process. 
The difficulty comes from the fact that the true CCQE 
events are masked by FSI in a detector built from nuclei. 
In general, at Cherenkov detectors such as MiniBooNE 



CCQE-like events are all those where no pion is detected 
while in tracking detectors such as K2K-SciBar/SciFi 
CCQE-like events are those where a single proton track 
is visible and at the same time no pions are detected. In 
both detector types the FSI lead to misidentified events, 
e.g., an initial A whose decay pion is absorbed or which 
undergoes "pion-less decay" can count as CCQE event 
— we call this type of background events "fake CCQE" 
events. We denote every event which looks like a CCQE 
event by "CCQE-like". 

To investigate the relation between the CCQE-like and 
true CCQE cross section, we show their ratio as a func- 
tion of proton and pion momentum thresholds in Fig. [3] 
As the proton is not at all relevant for the CCQE identi- 
fication in Cherenkov detectors, the ratio is independent 
of the proton momentum detection threshold (dashed line 
in left panel). This is very different in tracking detectors 
which rely on the detected proton — here the efficiency 
is reduced to «10% at a proton momentum threshold of 
0.5 GeV (soHd fine in left panel). Even at = 

the efficiency does not exceed 80% because of charge- 
exchange processes that lead to the emission of unde- 
tected neutrons and because of secondary proton knock- 
out that leads to multiple-proton tracks. Focussing on the 
right panel of Fig. [3] we find that the CCQE-like cross 
section increases for both detector types as |/?lthj.gs in- 
creases. In this case even more events with pions in the 
final state appear as CCQE-like because these pions are 
below threshold and thus not detected. 

Fixing the flux normalization with HARP's pion- 
production data, the MiniBooNE collaboration has pre- 
sented their first, preliminary absolutely normalized to- 
tal, differential and double differential cross sections for 
CCQE and finds an excess of about 35% compared to 
the total cross section measured at NOMAD, ANL and 
BNL 1 12]. We emphasize that these absolute cross sec- 
tions depend directly on the pion background subtraction 
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FIGURE 4. Double differential cross section on '^C averaged over the MiniBooNE flux as a function of the muon kinetic energy 
and the muon scattering angle. The left panel shows the true CCQE cross section, the right panel the ratio of the CCQE-like to the 
true CCQE cross section. 



whic h ag ain is based on the Monte Carlo prediction (cf. 
Ref. 01). 

In Fig. m we show our prediction for the double differ- 
ential cross section at MiniBooNE in muon observables, 
all calculated with = 1 GeV. The left panel shows 
the true CCQE events. To emphasize the role of "fake" 
CCQE events, we show the ratio CCQE-like/true CCQE 
in the right panel. This ratio shows that fakes contribute 
mainly at high energy transfers (low T^) and forward an- 
gles. Unlike for monochromatic beams, the QE and A 
peaks are not distinguishable any more but strongly over- 
lap. This fact makes a model-independent cut based on 
muon variables to subtract the background impossible. 
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MINIBOONE'S CC ln+ MEASUREMENT 

In Fig. |5]we give our results for the single- 7r+/QE ratio 
for CC interactions on mineral oil CH2. The solid lines 
denote the CCl;r+-Hke/CCQE-Hke result, the dashed 
lines stand for the true CCl;r+/true CCQE result, and 
the dash-dotted lines give the vacuum expectation. Note 
that we have applied the Cherenkov detector identifica- 
tion criteria. 

We emphasize that nuclear corrections cancel out in 
the ratio, only as long as FSI are not considered ("true" 
vs. "free"). In general, the complexity of FSI prevent 
such cancellations as one can infer from the result de- 
noted with "like" which does not coincide with the "true" 
and "free" ones. 

We further compare to very recent MiniBooNE data 
ifTsIl (Fig.|5]l. Let us first focus on the data denoted with 
the triangles (upper data set). These are corrected for FSI 
using a specific Monte Carlo generator, i.e., they give the 
cross sections for bound nucleons "before FSI". As this 
procedure introduces a model dependence in the data, a 



FIGURE 5. Single-;r+/QE cross section ratio for CC in- 
teractions vs. neutrino energy on CH2 together with recent 
data from MiniBooNE 1 13] (upper data set: corrected for FSI, 
lower data set: uncorrected for FSI). The solid lines denote 
the CCl;r+-like/CCQE-like result (Cherenkov detector defini- 
tions), the dashed lines stand for the true CCl;r+/true CCQE 
result, and the dash-dotted lines give the vacuum expectation, 
i.e., the sum of the nucleon cross sections (with two additional 
protons in the MiniBooNE case). 



fully consistent comparison is not possible. Ignoring this 
inconsistency, our calculation denoted with "true" should 
be the one to compare with. The agreement is perfect for 
energies up to 1.5 GeV, and still within their error bars 
for higher Ey ■ The MiniBooNE data denoted with bullets 
(lower data set) is their result for the ratio of CCl ;r+-like 
to CCQE-like. As these data are not corrected for FSI 
within a specific Monte Carlo event generator, this ob- 
servable is less model dependent. Still, the energy recon- 
struction requires specific assumptions as well as the de- 
tector simulation. We find that our calculation clearly un- 
derestimates the uncorrected data. However, the perfect 




FIGURE 6. Left panel: NC induced single-;r^ production on '^C as a function of tlie pion kinetic energy averaged over the 
MiniBooNE flux. Rigiit panel: same on '^O averaged over the K2K flux. The dashed lines show our calculation without FSI or 
spectral functions, both included in the full calculation denoted with the solid lines. The dash-dotted lines indicate the contribution 
from the A resonance to the full calculation. 



agreement with their corrected distribution indicates a 
significant difference between the pion absorption mod- 
els. 

The underestimate of the pion/quasielastic ratio in 
particular at higher energies could be due to, among other 
possibilities, an underestimate of the pion production 
cross section or an overestimate of the CCQE-like cross 
section. Both depend directly on the input at the nucleon 
level, i.e., in particular on the axial form factor C5 of the 
A resonance (see Fig. [Hi. For the results presented here, 
we have used the ANL data as a reference. However, we 
have shown in Ref. Iil4,l that an increase of the total pion 
production cross section on the nucleon compatible with 
the BNL data also seems to be insufficient to describe 
this ratio at all energies. We note that a similar result for 
the ratio has been recently obtained by Athar et al. ifTsll . 

NCl;r° 

In the left panel of Fig. |6] we show our results for NC 
single-;r'' production off '^C as a function of the pion ki- 
netic energy. We have averaged over the MiniBooNE en- 
ergy flux which peaks at about 0.7 GeV neutrino energy 
II16I1 . In NC reactions the total pion yield is dominated by 
Tt" production, while n'^ dominate in CC processes (for 
details, see Ref. iITtIi and references therein). Comparing 
the dashed with the solid line (results without FSI and 
spectral function vs. full calculation), one finds a consid- 
erable change. The shape is caused by the energy depen- 
dence of the pion absorption and rescattering cross sec- 
tions. Pions are mainly absorbed via the A resonance, i.e, 
through nN — > A followed by AA^ NN. This explains 
the reduction in the region around 0.1 — 0.3 GeV. 
Pion elastic scattering nN — > nN reshuffles the pions to 
lower momenta and leads also to charge exchange scat- 
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FIGURE 7. NC induced single-;r'' production on CH2 as a 
function of the pion momentum averaged over the MiniBooNE 
flux. The dashed and the solid lines denote the calculation with 
a modified dipole form factor (fitted to ANL data) for the A 
resonance (identical to the result in the left panel of Fig.[6l(. The 
dash-dotted and dotted lines are obtained with a dipole form for 
the form factor (fitted to the BNL data). The MiniBooNE data 
are taken from Ref. 11^ . 



tering into the charged pion channels. The vast majority 
of the pions comes from initial A excitation (dash-dotted 
line), their production in the rescattering of nucleons is 
not significant at these energies. 

The right panel of Fig. |6] shows the results for NC 
single- tt" production off ^^O averaged over the K2K en- 
ergy flux which peaks at about 1 .2 GeV neutrino energy 
II19I1 . Compared to the left panel, the spectrum is broader 
and extends to larger Tj^ due to the higher neutrino en- 
ergy. Again, pion production through initial QE scatter- 
ing is not sizable. 

In Fig. I2] we show also the latest NCtt" differential 



cross sections from MiniBooNE ifTsll . It can be seen that 
the full line, depicting the results of our calculation after 
FSI with the fit of the elementary cross section to the 
ANL data lies nearly a factor 2 below the data. The 
dashed line gives the result of the same calculation before 
FSI, the dash-dotted line shows the same quantity, but 
using a fit to the BNL data for the elementary cross 
section. We find it remarkable that these two curves 
(without FSI) are comparable to the data (with FSI), even 
though pion FSI have a major effect on the spectra (see 
Fig. |6ll. The absolute height of the experimental cross 
section is thus hard to reconcile with present descriptions 
of the elementary cross section and with what we know 
from other reactions about the pion FSI. 

DUALITY IN LEPTON SCATTERING ON 
THE NUCLEON 

At higher energies that correspond to invariant masses of 
the W,ZN system > 2 GeV the DIS channel becomes es- 
sential. A connection between the asymptotic DIS region 
and the resonance-dominated mass range is provided by 
the so-called Bloom-Gilman quark-hadron duality. 

Nearly forty years ago. Bloom and Oilman found 
||2Q1 that in electron scattering on protons the inclu- 
sive structure function F2 in the resonance region oscil- 
lates around the DIS scaling curve and, after averaging, 
closely resembles it. Furthermore, the resonance region 
data "slide" along the DIS curve with increasing Q^. If 
duality is understood quantitatively, there may be vari- 
ous applications. For example, the region of high Bjorken 
variable x is hardly experimentally investigated, because 
in the DIS region it would require very high and thus 
huge luminosities. If duality is satisfied with good accu- 
racy, one would be able to use the data in the resonance 
region to reach high x at reasonable 

The topic becomes even more interesting when turn- 
ing to nuclear targets and neutrino sources. The cur- 
rent precision measurements of the oscillation parame- 
ters require an efficient and accurate description of the 
neutrino-nucleus cross sections. Of particular interest is 
the resonance region and the possibility of linking it with 
the DIS region. A hadronic description of a neutrino- 
nucleus cross section at low Q~ requires a good knowl- 
edge of vector and axial transition form factors for each 
resonance. For the majority of the resonances, these tran- 
sition form factors are not well constrained. Provided that 
one can establish that quark-hadron duality holds with a 
reasonable accuracy, one could think of using the DIS 
results for estimating the neutrino-nucleus cross sections 
in the transition region. 

So far, most theoretical studies of quark-hadron du- 
ality in lepton scattering were dealing with nucleon tar- 



gets 112 IL I22I I23L 12411 . The DIS parts were considered as 
known, the structure functions in the scaling region be- 
ing conventionally evaluated from leading twist (LT) par- 
ton distribution functions (PDF): for example, F^'^^'"^^ = 
{F^'' + ')/2 = 5x/18 • (m + M + + J + 2^/5 + 2s/5), 
pVN{LT) ^ ^pvp ^ ^v„) i2=x{u + u + d + d + s + s). 
For nucleons, several parametrizations of the PDFs are 
generally available (from the ORV, CTEQ and MRST 
groups). In the region of moderate x, which is of inter- 
est for our duality study, they provide nearly the same 
results. 

The studies of the resonance region differ in the way 
the models treat the resonant contributions and the way 
they extract the structure functions. 

The fact, that the resonance data for increasing Q2 
"slide" along the DIS curve, can be observed if a few |24] 
or even only one ll2lll resonance are taken into account. 
The advantage of the model fl^ is that the structure 
functions are given as simple analytical functions of the 
momentum transfer squared and the energy transfer 
V, provided that the resonance form factors are known. 
In this work the first four resonances were considered. 
Oenerally, however, as it was argued by Close f^S*], 
inclusion of several resonances of different parities is 
desirable. 

A quark model for resonance excitation has been ap- 
plied by |22, 23] for the investigation of the duality in 
electro- and neutrinoproduction. 

Within the OiBUU framework, as it was mentioned 
above, 13 resonances can be considered for both electron 
and neutrino reactions. The cross section is calculated 
numerically and the structure function F2 is extracted 
from the cross section in a convenient way: 
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taken from ll26ll . The coefficients 
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are the Mott cross sections for electron and neutrino 
reactions. For electroproduction on an isoscalar target, 
C7^ = {a'^'^ + (7™)/2 is half sum of electroproduction 



cross sections on proton and neutron. For charged cur- 
rent neutrinoproduction, in order to eliminate the struc- 
ture function Ft,, one should use the linear combination 
of the neutrino and anti-neutrino cross sections. For an 
isoscalar target it is sufficient to take C7^ = {g^'' + C7^'' + 

C7^« + C7^'')/4. 

For a quantitative estimate of the validity of duality it 
is convenient to introduce the ratio of the integrals of the 
resonance (res) and DIS structure functions 



2 ] 
DIS' 
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where denotes 2xFi , F2 or xFj, (for neutrino scatter- 
ing). The value 2d/s is taken as the actual value for a 
given parametrization of DIS PDFs (for nucleon, in our 
case, = 10 GeV) or a DIS experimental data set (for 
nuclei). Under conditions of perfect quark-hadron dual- 
ity this ratio would be 1 and independent of Q^. Thus, the 
degree to which the local duality is fulfilled can be esti- 
mated from the dependence and the deviation from 1 
of the computed 

The isoscalar nucleon structure function, which 
includes both resonance and background contributions, 
is shown in Fig. [8] versus the Nachtmann variable ^. 
Notice, that ^ decreases with increasing invariant mass 
W. For a given value, the highest peak at the larger 
£, value corresponds to the A— resonance peak, and the 
two lower peaks at smaller values of ^ correspond to 
the second (1.40 GeV < W < 1.56 GeV) and the third 
(1 .56 GeV < W < 2.0 GeV) resonance regions. The gen- 
eral picture shows a reasonable agreement with the dual- 
ity hypothesis. 
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FIGURE 8. as a function of t, , for = 0.225, 0.525, 

1.025 and 2.025 GeV^ (indicated on the spectra), compared 
with the leading twist parametrizations at Q = 10 GeV^. 

In Fig. |9] the ratio of the integrals if^, defined in (O, 
is shown not only for the whole structure function (reso- 
nance + 1-pion background), but also for the resonance 
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region to the leading twist functions. 



contribution separately. For > 0.5 GeV^, the ratio if 
for the resonance contribution only is at the level of 0.85, 
which is smaller and flatter in in comparison with the 
results I24ll27il of the Dortmund group resonance model. 
The difference is due to the different parametrization of 
the electromagnetic resonance form factors used in the 
two models. The In background gives a noticeable con- 
tribution and brings the ratio up to 0.95. The fact, that 
the ratio is smaller than 1 is no surprise, because addi- 
tional non-resonant contributions like 2- and many-pion 
background are possible, but not taken into account here. 
They are the subject of coming investigations. Indeed, 
they could be obtained from requirement /2 = 1 ■ 

The principal feature of neutrino reactions, stemming 
from fundamental isospin arguments, is that duality does 
not hold for proton and neutron targets separately. The 
interplay between the resonances of different isospins al- 
lows for duality to hold with reasonable accuracy for the 
average over the proton and neutron targets. We expect, 
that a similar picture emerges in neutrino reactions with 
nuclei. 

For neutrinoproduction, the structure function fj^^ is 
shown in Fig.[TO]for the resonance contribution only. 

The ratio /j ^ is shown in Fig. [TT] and appears to be 
at the level of 0.7, which is (similar to the electron case) 
smaller than 0.8, which has been calculated within the 
Dortmund resonance model ll24l IZT^. Thus, one would 
expect a large contribution from the background. The 
role of the background in the anti-neutrino channel is 
under investigation. 




FIGURE 10. F^^ as a function of t, , for = 0.225, 0.525, 
1.025 and 2.025 GeV^ (indicated on the spectra), compared 
with the leading twist parametrizations at 2 = 10 GeV^. 
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FIGURE 11. Ratio i:^^ of the integrated F^^ in the reso- 
nance region to the leading twist functions. 



GiBUU model. The FSI do not play any role for the in- 
clusive cross section. The in-medium corrections are of 
considerable importance. 

In Fig. [12] the resonance contribution to the F2 /A 
structure functions for a carbon target are shown for sev- 
eral values. They are compared to experimental data 
obtained by the BCDMS collaboration lUflHI)] in muon- 
carbon scattering in the DIS region (2^ ~ 30 — 50 GeV^). 
They are shown as experimental points connected by 
smooth curves. For different values, the experimen- 
tal curves agree within 5% in most of the t, region, as 
expected from Bjorken scaling. 

When investigating duality for a free nucleon, we took 
the average over free proton and neutron targets, thus 
considering the isoscalar structure function. Since the 
carbon nucleus contains an equal number of protons and 
neutrons, averaging over isospin is performed automati- 
cally. Due to the Fermi motion of the target nucleons, the 
peaks from the various resonance regions, which were 
clearly seen for the nucleon target, are hardly distinguish- 
able for the carbon nucleus. 
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DUALITY IN LEPTON SCATTERING ON 
THE NUCLEUS 

Recent electron scattering measurements at JLab have 
confirmed the validity of the Bloom-Gilman dual- 
ity for proton, deuterium and iron structure functions. 
Further experimental efforts are required for neutrino 
scatterin g. A mong the upcoming neutrino experiments, 
Minerva ll28ll and SciBooNE[29] aim at measurements 
with carbon, iron and lead nuclei as targets. 

One of the major issues for nuclear targets is the def- 
inition of the nuclear structure functions F^^^^y Exper- 
imentally they are determined from the corresponding 
cross sections, using Eq. (|2|l. 

We follow the same procedure, using the GiBUU cross 
sections. So, at the first step the inclusive double differ- 
ential cross section da jdQ^dv is calculated within the 



FIGURE 12. Resonance curves F2 ''/12 as a function of ^, 
for = 0.45, 0.85 and 1.4 GeV^ (indicated on the spectra), 
compared with the experimental BCDMS data L30i . i31i1 in the 
DIS region at 2^^^ = 30, 45 and 50 GeV^. 

As expected from local duality, the resonance structure 
functions for the various values slide along a curve, 
whose ^ dependence is very similar to the scaling-limit 
DIS curve. However, for all ^, the resonance curves lie 
below the experimental DIS data. 

To quantify this underestimation, we now consider the 
ratio of the integrals of the resonance (res) and DIS struc- 
ture functions, determined in Eq. ([3]). As it is explained in 
il^, the integration limits are to be determined in terms 
of the effective W variable, experimentally (see, for ex- 
ample, |33]) defined as — + ImnV — Q^. For a 
free nucleon W coincides with the invariant mass W. For 
a nucleus, it differs from W due to the Fermi motion of 
bound nucleons, but still gives a reasonable estimation 



for the invariant mass region involved in the problem. 

In particular, the resonance curves presented in all fig- 
ures are plotted in the region from the pion-production 
threshold up to W = 2 GeV. For a free nucleon, the 
threshold value for 1-pion production (and thus the 
threshold value of the resonance region) is Wmin = 
Winin ~ 1 . 1 GeV. Bound backward-moving nucleons in 
a nucleus allow lower W values beyond the free-nucleon 
limits. The threshold for the structure functions is now 
defined in terms of v or W, rather than W. Hence, we 
consider two different cases in choosing the ^ integration 
limits for the ratio (|3]l. First, for a given Q^, we choose 
the £, limits in the same manner as for a free nucleon: 

^mm = <^(W = 1.6GeV,e2)^ 

(4) 

^max = ^(W = l.lGeV,e2) . 

We refer to this choice as integrating "from 1.1 GeV". 
The integration limits for the DIS curve always corre- 
spond to this choice. As a second choice, for each 
we integrate the resonance curve from the threshold, 
that is from as low W as is achievable for the nucleus 
under consideration. This corresponds to the threshold 
value at higher £, and is referred to as integrating "from 
threshold". With this choice we guarantee that the ex- 
tended kinematical regions typical for resonance produc- 
tion from nuclei are taken into account. Since there is no 
natural threshold for the ^,„,„, for both choices it is deter- 
mined from W = 1.6 GeV, as defined in Eq. 

The results for the ratio ([3j are shown in Fig.[T3] The 
curve for the isoscalar free-nucleon case is the same as 
in Ref. |24] with the "GRV" parametrization for the DIS 
structure function. One can see that the carbon curve ob- 
tained by integrating "from threshold" lies above the one 
obtained by integrating "from 1.1 GeV", the difference 
increasing with Q^. This indicates that the threshold re- 
gion becomes more and more significant, as one can see 
from Fig. [13] Recall, that the flatter the curve is and the 
closer it gets to 1, the better local duality would hold. Our 
calculations for carbon show that the ratio is at the same 
level as that for the free nucleon or even higher. 

For neutrino-iron scattering, the structure functions 
F2^'^ are shown in Fig. [14] As for the electron-carbon 
results of Fig. [12] the resonance structure is hardly vis- 
ible. The resonance structure functions are compared to 
the experimental data in the DIS region obtained by the 
CCFR |34] and NuTeV |35] coflaborations. It appears, 
that the resonance curves slide along the DIS curve, as 
one would expect from local duality, but lie well below 
the DIS measurements. Hence, the computed structure 
functions do not average to the DIS curve. The necessary 
condition for local duality to hold is thus not fulfilled. 

As expected from local duality, the resonance structure 
functions for the various values slide along a curve, 
whose £, dependence is very similar to the scaling-limit 
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FIGURE 13. Ratio defined in Eq.l[3]l for the free nucleon 
(dash-dotted line), and '^C. We consider the lower limits deter- 
mined by ly = 1.1 GeV (solid line) and by the threshold value 
(dotted line). 
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FIGURE 14. F^^'^yse as a function of ^ for = 0.2, 0.45, 
0.85, and 1.4 GeV^. The calculations are compared with the 
DIS data from CCFR and NuTeV experiments Refs. Isllsll. 
The DIS data refer to measurements at = 7.94, 12.6 and 
19.95 GeV^. 

DIS curve. However, for all the resonance curves lie 
below the experimental DIS data. 

The ratio /j defined in Eq.© is shown in Fig. [15] 
The curve for the isoscalar free nucleon case is also pre- 
sented for comparison. Our results show, that 1 ) this ratio 
is significantly smaller than 1 for all Q^; 2) it is signifi- 
cantly smaller than the one for the free nucleon; 3) I2 is 
even lower than the corresponding ratio for electropro- 
duction; 4) I2 slightly decreases with Q^. 

To summarize, we find that the resonance structure 
functions are consistently smaller than DIS functions in 
the same region of the Nachtmann variable ^. This is 
in agreement with earlier work [32], which implements 
elementary resonance vertices and nuclear effects differ- 
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FIGURE 15. Ratio ^ defined in Eq. lO for the free 
nucleon (dasii-dotted line) and ^^Fe. The results are displayed 
for two choices of the lower limit in the integral: W = l.l GeV 
(solid line) and threshold (dotted line). For each of these two 
choices we have used two sets of DIS data in determining the 
denominator of Eq. (O. These sets of DIS data are obtained at 
qIjs = 12.59 and 19.95 GeV^. 



ently. Recall that in this analysis for nuclei we include the 
resonance structure functions and ignore the background 
ones. To estimate their contribution and compare the re- 
sults with the nucleon case is one of the primary tasks of 
coming investigations. 



CONCLUSIONS 

Overall, the impact of nuclear effects impact on observ- 
ables is dramatic. Final state interactions of produced pi- 
ons lead to pions being absorbed in the nucleus or rescat- 
tered, thus reducing the cross section and shifting the ob- 
served pion distributions to lower values of the pion en- 
ergy. Furthermore, such interactions can lead to knock- 
out nucleons which experiments tend to identify with 
QE events. CCQE and pion production are thus closely 
entangled. Thus model-independent data are needed for 
a meaningful comparison of theory with current exper- 
iments. Initial state interactions reduce the cross sec- 
tions and structure functions in the resonance region 
and change the form of the distributions, dissolving the 
peaks of individual resonances. The Bloom-Gilman du- 
ality seems to be violated for nuclear targets, possibly 
indicating a major change of background amplitudes in 
the medium. 



ACKNOWLEDGMENTS 

This work has been supported by the Deutsche 
Forschungsgemeinschaft. 



REFERENCES 

1. T. Leitner, O. Buss, L. Alvarez-Ruso and U. Mosel, Phys. 
Rev. C79, 034601 (2009), arXiv:0812.05^ 

2. MAID, http : / /www . kph . uni-mainz . de/ MAID[ 

3. L. Tiator and S. Kamalov, arXiv:nucl-th/0603012 

4. D. Drechsel, S. S. Kamalov and L. Tiator, Eur. Phys. J. 
A34, 69 (2007), arXiv:07 10.0306 

5. S. J. Barish et al, Phys. Rev D19, 2521 (1979). 

6. G. M. Radecky etai, Phys. Rev. D25, 1161 (1982). 

7. T. Kitagaki et al, Phys. Rev. D34, 2554 (1986). 

8. GiBUU, 

jhttp : / /gibuu ■ physik .uni-giessen. de/GiBUU[ 

9. B. Krusche et al, Eur. Phys. J. A22, 277 (2004), 
arXiv:nucl-ex/0406002 

10. O. Buss, L. Alvarez-Ruso, P. Muehlich and U. Mosel, 
Eur. Phys. J. A29 (2), 189 (2006), arXiv:nucl-th/0603003] 

11. O. Buss, L. Alvarez-Ruso, A. B. Larionov and U. Mosel, 
Phys. Rev. C 74, 044610 (2006), arXiv:nucl-th/0607016j 

12. MiniBooNE, T. Katori, arXiv:0909.1996 [hep-ex]. 

13. MiniBooNE, A. A. Aguilar-Arevalo et al., Phys. Rev. 
Lett. 103, 081801 (2009), arXiv:0904.3159 

14. T. Leitner, O. Buss, U. Mosel and L. Alvarez-Ruso, 
arXiv:0909.0838 

15. M. S. Athar, S. Chauhan and S. K. Singh, 
arXiv:0908. 14421 

16. MiniBooNE, A. A. Aguilar-Arevalo et al., Phys. Rev. 
D79, 072002 (2009), arXiv:0806.1449 

17. T. Leitner, O. Buss, U. Mosel and L. Alvarez-Ruso, Phys. 
Rev. C79, 038501 (2009), arXiv:0812.1787 

18. C. Anderson, Neutral current In^ production 
at MiniBooNE, talk given at NUINT 09, May 
17-22, 2009, Sitges, Spain, available online at 
http://nuint09.ifae. es /Agenda . html 

19. K2K, S. Nakayama et al., Phys. Lett. B619, 255 (2005), 
arXiv:hep-ex/0408134 

20. E. D. Bloom and F. J. Oilman, Phys. Rev. Lett. 25, 1 140 
(1970). 

21. K. Matsui, T. Sato and T. S. H. Lee, Phys. Rev. C72, 
025204 (2005), arXiv:nucl-th/0504051 

22. K. M. Graczyk, C. Juszczak and J. T. Sobczyk, Nucl. 
Phys. A781, 227 (2007), arXiv:hep-ph/0512015 

23. K. M. Graczyk, arXiv:0910.2355 [hep-ph]. 

24. O. Lalakulich, W. Melnitchouk and E. A. Paschos, Phys. 
Rev. C75, 015202 (2007), arXiv:hep-ph/0608058 

25. R E. Close and W. Melnitchouk, Phys. Rev. C68, 035210 
(2003), arXiv:hep-ph/0302013 

26. L. W. Whiflow, E. M. Riordan, S. Dasu, S. Rock and 
A. Bodek, Phys. Lett. B282, 475 (1992). 

27. O. Lalak ulich etal, AIP Conf. Proc. 967, 24 3 (2007). 

28. Minerva, ll http : //minerva . f nal . gov| 

29. SciBooNE, http : / /www- sciboone . ifna l . gov/l 

30. D. Bolhni et al., Phys. Lett. B104, 403 (1981). 

31. BCDMS, A. C. Benvenuti et al, Phys. Lett. B195, 91 
(1987). 

32. O. Lalakulich, N. Jachowicz, C. Praet and J. Ryckebusch, 
Phys. Rev C79, 015206 (2009), arXiv:0808.0085 

33. R. M. Sealock et al., Phys. Rev. Lett. 62, 1350 (1989). 

34. W. G. Sehgman etal, Phys. Rev Lett. 79, 1213 (1997). 

35. NuTeV, M. Tzanov etal, Phys. Rev D74, 012008 (2006), 
|arXiv:hep-ex/0509010, 



